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Abstract 

The strong convergence of Euler approximations of stochastic delay differential equations is 
proved under general conditions. The assumptions on drift and diffusion coefficients have been 
relaxed to include polynomial growth and only continuity in the arguments corresponding to 
delays. Furthermore, the rate of convergence is obtained under one-sided and polynomial Lip- 
schitz conditions. Finally, our findings are demonstrated with the help of numerical simulations. 

Keywords: stochastic delay differential equations, Euler approximations, strong convergence, 
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1 Introduction 

In modeling many real world phenomena, the future states of the system depend not only on 
the present state but also on its past state(s). The models based on stochastic delay differential 
equations (SDDEs) could be used in such situations. They have found numerous applications in 
various fields, for example, in communications, physics, biology, ecology, economics and finance. 
One could refer to pQ - |1U] and references therein. 
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In many applications, one finds SDDEs for which solutions can not be obtained explicitly and 
hence this necessitates the development of numerical schemes. Over the past few years, many 
authors have shown their interest in studying strong convergence properties of numerical schemes 
of SDDEs, to name a few, [11] - [16j . Recently there has been a growing interest in weak convergence 
of numerical schemes as well, see for example, [T7], [18] and |19] , We also mention a paper by |20] 
on pathwise approximations and its relation to strong convergence. For the purpose of this article, 
we study strong convergence of Euler approximations of SDDEs due to its wide applicability. We 
exploit the fact that SDDEs can be regarded as special cases of stochastic differential equations 
(SDEs) with random coefficients. To this end, first we establish strong convergence results of 
Euler schemes of SDEs with random coefficients under mild conditions. This, in turn, allows us to 
prove the convergence in the case of SDDEs under more relaxed conditions than those existing in 
the literature. More precisely, to the best of our knowledge, strong convergence has been proved 
so far by assuming linear growth and local Lipschitz conditions in both variables corresponding 
to delay and non-delay arguments, see for example [16J. Whereas, we relax these conditions by 
assuming polynomial growth and continuity in the argument corresponding to delays. Moreover 
no smoothness condition on the initial data is assumed. Therefore, the setting in this article is 
more general than those present in the literature, see for example, |llj . |14j . [15j . [16] and [20] , In 
addition, it is shown here that the rate of convergence is one-fourth when the drift coefficient satisfies 
one-sided Lipschitz and polynomial Lipschitz in the non-delay and delay arguments respectively. 
This result is in agreement with the findings of [21j . 

We conclude this section with the introduction of some basic notation. The norm of a vector 
x G M rf and the Hilbert-Schmidt norm of a matrix A 6 M rfxm are respectively denoted by \x\ and 
\A\. The transpose of a vector x S M d is denoted by x T and the inner product of two vectors 
x,y G M. d is denoted by (x,y) := x T y. The integer part of a real number x is denoted by \x\. We 
denote the indicator function of set A by 1a- Moreover, C p = C p {yt,F, P) denotes the space of 
random variables X with a norm \\X\\ p := (E[|X| P ]) 1 ^ P < oo for p > 0. Finally, V denotes the 
predictable <r-algebra on M + x and B(V), the cr-algebra of Borel sets of topological spaces V. 
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2 Main Result 



Let (O, {J~t} {t>o} > J~ ,¥) be a filtered probability space satisfying the usual conditions, i.e. the filtra- 
tion is increasing and right continuous. Let {W (t)} {t>o} be an m-dimensional Wiener martingale. 
Furthermore, it is assumed that (3(t,yi, ... ,yf,,x) and a(t,y\, ... ,yk,x) are <8> B(R. dxk ) ® 

£>(R d )-measurable functions and take values in M rf and M cixm respectively. For a fixed T > 0, let 
the stochastic delay differential equation (SDDE) on (0, {J^jj^o}, J 7 , P) be defined as follows, 

dX{t) = P{t,Y(t),X(t))dt + a(t,Y{t),X(t))dW{t), te[0,T], 
X(t) = £(t), *G[-^0], 

where {f(t) : -# < t < 0} G C^ ([-fl", 0]; R d ) for some if > and Y(t) := (X((5i(*)), . . . , X((5 fc (t))). 
The delay parameters <5i(t), . . . ,5f;(t) are increasing functions of t and satisfy —H < 5j(t) < [^]r 
for some r > and j = 1, . . . , k. 

Further, for every n > 1 and t G [0, T], the Euler scheme of SDDE (pQ) is given by 



(2) 



dX n (t) = (3(t,Y n (t),X n (K n (t)))dt + a(t,Y n (t),X n (K n (t)))dW(t), t e [Q,T], 
X n (t) = e(t), «G[-fT,0], 

where ^(t) := (X n (6i(t)), . . . , X n (5k(t))) and n n is defined by 

n n {t),= ^tlM +t0 (3 ) 
n 

with the observation that for equation ([2]) one takes to = in ([3]). We note that two popular cases 
of delay viz. 5i(t) = t — r and <5j(t) = | r can be addressed by our findings. This type of delay 
parameters have been extensively discussed and found wide applications in the literature. The 
reader could consult the following, for example, pQ - [6], [10] and references therein. 

Remark 1. We note that the Euler scheme ([2]) defines approximations to SDDEs in an explicit way 
without a discretization of the delay terms. From our main theorem and corollaries on convergence 
of scheme ([2]), one could easily obtain results on convergence of this scheme with discretized delay 
terms. Therefore for matters of notational simplicity, we choose the former approach. 

Let y := (yi, . . . ,yk). We make the following assumptions for our result. 



C-l. There exist constants G > and I > such that for any t G [0, T], 



\P(t,y,x)\ + \a{t,y,x)\ < G(l + \y\ l + \x\) 

for all x eR d and y e R dxk . 

C-2. For every R > 0, there exists a constant Kr > such that for any t G [0,T], 

(x - z, p(t, y, x) - pit, y, z)) V \a(t, y, x) - a(t, y, z)\ 2 < K R \x - z\ 2 
whenever \y\, \z\ < R. 

C-3. The functions /3(t, y, x) and a(t, y, x) are continuous in y uniformly in x from compacts, i.e. 
for every R > and t G [0, T], 

sup{\P(t,y,x)-P(t,y',x)\ + \a(t,y,x)-a(t,y',x)\}->0 as y ->• y. 

\x\<R 

The conditions C-l and C-2 are sufficient for existence and uniqueness of solution of SDDE ([T]) 
and the Euler scheme ([2]) (see [12], [22]). We state the main result of this paper in the following 
theorem. 

Theorem 1. Suppose C-l to C-3 hold, then the Euler scheme ([2|) converges to the true solution of 
SDDE © in C p -sense, i.e. 



lim E 



sup \X(t)-X n (t)\P 

0<t<T 



for all p > 0. 



Remark 2. We note that we can assume, without loss of generality, that T is a multiple of r. If 
not, then SDDE ([1]) and its EM scheme ([2]) are defined for T' > T with T' = Nt, where N is a 
positive integer. The results given in this report are then recovered for the original SDDE ([1]) by 
choosing drift and diffusion parameters as P1u<t} an d a ^{t<T}- 

First, we develop requisite theory for SDEs with random coefficients in the following section 
and then prove Theorem Q] in Section [H 
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3 SDEs with Random Coefficients 

Let b(t, x) and er(t, x) be V (g> Z^R^-measurable functions which take values in M. d and M dxrn 
respectively. Suppose < to < ti < T, then let us consider an SDE with random coefficients given 
by 

dX(t) = b(t, X(t))dt + <x(i, X(t))dW(t), V t e [to,h], (4) 

with initial value X(to) which is an almost surely finite J-£ -measurable random variable. 

For every n > 1, let b n (t,x) and a n (t,x) be V ® i3(M (i )-measurable functions which take values 
in W d and ^ dxm respectively. Then for any t G [to,ti], we define the following Euler scheme 
corresponding to the SDE @ as, 

dX n {t) = b n {t, X n {K n (t)))dt + a n (t, X n {n n {t)))dW{t) (5) 

with initial value X n (to), which is an Tt -measurable, almost surely finite, random variable and n n 
as defined in ([3|). 

We make the following assumptions. 

A-l. There exists a constant c > and non-negative random variables {M n } n >i and M with 
E[M'P] V sup n >! E[M%] < N for every p > and some N := N(p) > such that, almost surely, 

|&(*,s)| + W(t,x)\ < c(M + \x\) 
\b n (t,x)\ + \a n (t,x)\ < c(M n + |x|) 

for any t G [0, T] and x G R d . 

A-2. For every R > 0, there exists an Tt Q -measurable random variable Cr such that, almost surely, 
for any t G [to, tl], 

(x - y, b{t, x) - b(t, y)) V \a(t, x) - a(t, y)\ 2 < C R \x - y\ 2 
for all |y| < R. Moreover, there exists a function / : R + — > M + such that 
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lim ¥(C R > f(R)) = 0. 

R— >oo 



A-3. For every R > and p > 



lim E 



Jt, 



t \x\<R 



sup {\b n (t,x) - b(t,x)\ p + \a n (t,x) - a(t,x)\ p }dt = 0. 



A-4. The initial values of SDE Q and its Euler scheme ([5]) satisfy 



lim E[\X(t )-X n (t )\P] =0 



for every p > 0. 

A-5. The initial values of SDE @ and its Euler scheme ([5]) have bounded p-th moments, i.e. there 
exists L := L(p) > such that 



for every p > 0. 

We note that assumptions A-l and A-2 are sufficient for existence and uniqueness of solution 
of SDE (HD (see [23J and [24|). Before proving the main result (Theorem [2]) of this section, we 
establish some lemmas. 

Lemma 1. Suppose that A-l and A-5 hold, then for some K := K(p,T,L) > 0, 



n\X(t )\ p ] V supE[|X n (t )H < L 



n>l 



E sup \X(t)\ p < K 



t <t<t x 



and 



supE sup \X n (t)\ p < K 



(6) 



n>l L f <*<*i 



for every p > 0. 



Proof. First one chooses any p > 2 



Then for the true solution {X(t)} t€ [ to tl j, by using Holder's 
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inequality, one obtains, 



\X(t)\P < 3P- 1 i|X(t )| p + |t-t r 1 [ t \b(s,X(sWds+ f 

y J to J to 



a(s,X(s))dW(s) 



Thus on taking supremum over [to, u] for some u G [to, h] and by applying Burkholder-Davis-Gundy 
inequality with constant c := c(p) > 0, one observes that 



E 



sup \X(t)\ p 

to<t<U 



< 3 p - 1 |E[|X(t )n +T p - 1 E 



|6(s,X(s))| p ds 



to 



+ cE 



|<r(s,X(s))| p ds 



L Jt 



Now A-l and A- 5 give 



E 



sup |X(t)| p 

to<t<U 



< 2 p - 1 |E[|X(t )| p ] + [T^ + cpE ^ n {|M| + |X( S )|} P (i S 

< 3 p - 1 |E[|X(t )| p ] + 2 p " 1 [r p - 1 + cjc^TA^ 



J to 



sup |X(r)| p 



which on the application of Gronwall's inequality yields 



E 



sup \X(t)\ p 

t <t<tl 



For any p > 2, we adopt similar arguments for {A^ n (t)}{t g [i 0<1 ]} in order to prove (jfij) with the 
remark that A- 5 guarantees that K does not depend on re. Then one observes that the application 
of Holder's inequality completes the proof. □ 

For every R > and re > 1, let us define the stopping times, 
t r := mf{t > t : \X(t)\ > R}, a nR := M{t > t : |X n (t)| > R} and v nR := r R A a nR , (7) 



where inf = oo. 

We prove below a very useful lemma of this article as it is used in the penultimate section in 
order to recover the rate of convergence. 
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Lemma 2. Let us assume that A-l and A-5 hold. Then 



lim E 

n— >oo 



\X n (s A v nR ) - X n (K n {s A v nR ))\ p ds 



to 







for every p > 0. 

Proof. In order to prove the result, one chooses p > 2. Then one immediately writes, 



E 



|X n (s A i/ nfl ) - X n (K n (s A twj))| p <is 



to 









= E 


[/" 


j 






J K 



"sAu nR 



' K n (sAv nR ) 



b n (r,X n (K n {r)))dr 



"SAv n R 

+ I a n (r,X n (K n (r)))dW(r) 

I K„(sAi/„ H ) 



ds 



which implies on applying Holder's inequality 



E 



ti 



\X n (s A v nR ) - X n (K n (s A u nR ))\ p ds 
ti 



to 

< 2 p_1 E 



sAv nR 



s /\ v n R - Kn(s A is nR )\ p M / | b n (r, X n (« n (r) )) \ p dr ) ds 

to \JK n (sAu nR ) 
rti 

-V 1 / E 

Jto 



sAu nR 



■ (-) / E 



-2 P_1 / E 
Jto 



a n (r,X n (K n (r)))dW(r] 

K n (sAv nR ) 

sAu nR 

\b n (r,X n ( Kn (r)))\ p dr ds 

K n (sAv nR ) 
ti r r sAu nR 



ds 



a n (r,X n ( Kn (r)))dW(r] 



' K n {sAv nR ) 

One finds bounds for the integrand of the first term of ([8]) as follows, 



ds. 



E 



sAu„ R 
{sAu nR ) 

< c p 2 p - 1 E 

< c p 2 p - 1 E 



|6 n (r,X n ( Kn (r)))|^r 



< c^E 



(8) 



Kn(sA^„H) 



(|M n | + |X n ( Kn (r))|fdr 



sAu nR 



(\M n \ p + \X n (K n (r))\ p )dr 



K, n {sAu nR ) 

\M n \ p + sup \X n (r)\ p )\s f\v nR -K n {s f\v nR ) 

to<r<t! 



< 



< 



n 

c p 2 p - x 



sup \X n (r)\ p 
t <r<ti 
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E[|M n |PJ +E 
(K + N). (using Lemma d]) 



(9) 
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For the second term of (JSj) , one writes, 



E 



sAu nR 



Kn(sAv nR ) 



a n {r,X n (n n {r)))dW{r) 



E 



Kn(s) 



1 [Kn (sAu nR ),sAu nR ] On (r, X n (k u {r)))dW{r) 



< i P -^lY /2 \s- Kn (s)\^E 



K n (s) 



1 [«n(sA 1 , nfl ), S A 1 , n fl]kn(r,X ri (K ri (r)))| p (ir 



which on the application of A-l and Lemma [T] yields 



E 



sAu nR 



K n (sAv nR ) 



a n (r,X n (K n (r)))dW(r) 



< I P(P -1^(?^-\K + N). (10) 



Substituting © and ^TUJ) in (jSJ yields 



E 



\X n {s A u nR ) - X n (K n (s A v nR ))\ p ds 



T(K + N)4P- l c p f 1 | (pip -I) 



n p/2 

< 0(n- p / 2 ). 



p/2 



p/2 



(11) 



Then by letting n — > oo on both sides, one obtains 



lim E 

n— ¥oo 



\X n {s f\ v nR ) - X n {n n {s A v nR ))\ p ds 



to 



for all p > 2. Then one observes that the application of Holder's and Jensen's inequalities completes 
the proof. □ 



In order to prove the following theorem, a similar approach to the one adopted by [25j for SDEs 
is followed here, but in a more general context. 

Theorem 2. Let A-l to A-5 hold and suppose there exists a unique solution {X n (i)}r to < t<tl } of 
the Euler scheme ([5]). Then, the Euler scheme © converges to the SDE @ in C p -sense, i.e. 



lim E 



sup \X(t)-X n (t)\ p 
to<t<h 



(12) 



for all p > 0. 



Proof. For every R > and n > 1, we consider stopping times as defined in ([7]). First we fix p > 4. 
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Then 



E 



sup \X(t) - X n (t)\P 

t <t<tl 

+E 



< E 



sup \X(t) - X n (t)\ p t {TR < tl or < tl or c R >f(R)} 
< <t<ti 



sup \X(t A u nR ) - X n (t A v n R)\ p ~k{C R <f(R)} 
t <t<tl 



(13) 



To estimate the first term of (fl3l) . one can use Young's inequality for q > p (1/p+l/g = 1) and 
T) > to obtain the following. 



E 



SUp \X(t) - X n (t)\ P l{ TR < tl or ffnK < tl or C R >f(R)} 
t <t<ti 



<^E 



sup \X(t)-X n (t)\* 

< <t<tl 



+- 



< — 2 9_1 < E 



sup 

t <t<*i 



+ E 
9-P 



3>/(q-p) 



< h oia nR < ti or C R > f(R)) 



sup |X n (t)| 9 

*0<i<*i 



qr) 



ipI (q-p) 



„ „,_ p) {p(tr < ti) + PKij < h) + ¥(C R > f(R))] 

+ ¥(C R > f(R))\, 



|e 


~\X(t r )\v~ 


+ E 


'\X n {a nR )f 




Rp 


RP 



which becomes 



E 



SUp \X(t) - X n (t)\ p t{ TR < tl or a nR <ti or C R >f(R)} 



<™2m + ^^\ 2 -^ + ncR>f{R)) 



•p/{i-p) I RP 



(14) 



Further for estimating the second term of (fT3l) . let e n (s) := X(s) — X n (s) for s G [toi*l]> then one 
could write 

de n (s) = b n (s)ds + a n (s)dW(s) 

where b n (s) := b(s, X(s)) - b n (s, X n (K n (s))) and a n (s) := a(s, X(s)) - a n (s, X n (n n (s))). Thus, on 
the application of Ito's formula, one obtains 

m m 

d\e n {s)\ 2 = 2(e n (s)X(s))ds + Y,K(s)\ 2 ds + 2^(e n (s),a i n (s))dW\ S ) 

-i=i i=i 
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which on integrating over s G [to,t A v n R\ for some t € [to,ii] yields 



m ftAu nR 



|e n (tA^ nJ? )| 2 = | en (t )| 2 + 2/ (e n (s), 6 n (s))ds + J2 K( s )f ds 

™ rtAu nR 
+2^ / (e n ( S ),4( S ))d^( S ) 
i=i Jt ° 

where a % n and W l denote the i-th column of (d x m)-matrix a and the i-th element of Wiener 
(to x l)-vector W respectively for all i £ {l,...,m}. Further, after raising the power to p/2 
for some p > 4 and applying Holder's and Burkholder-Davis-Gundy inequalities, one obtains the 
following, 



E 



SUp \e n (t /WnR)\ P t{C R <f(R)} 

tg<t<U 



< 2 p ~ 2 E[\e n {t c 



+2(8T) E 2 1 E 



1 {C R <f(R)} I \(e n (s),b n (s))\2ds 
t 



+(4mT)V J^E 



i=l 



uAu„ R 

Mc^KR)} I K(s)\ P ds 
to 



+2— m— T— c^E 



uhv nR 

1 {Cn<f(R)} I \(e n (s),a l n (s))\2ds 
to 



= 1 + 11 + 111 + IV 



(15) 



where c := c(p) is the constant of the Burkholder-Davis-Gundy inequality. 

In order to estimate 77, one observes that the application of A-2 and Cauchy-Schwarz inequality 
yields, 



(s)) = (X(s) - X n (K n (s)),b(s,X(s)) - b(s,X n (K n (s)))) 

+{X(s) - X n (K n (s)),b(s,X n (K n (s))) -b n (s,X n (K n (s)))) 

+{X n ( Kn {s)) - X n (s), b(s, X(s)) - b n (s, X n (K n (s)))) (16) 
< C R \X(s) -X n ( Kn (s))\ 2 + \X(s) - X n (K n (s))\\b(s,X n (K n (s))) - b n (s,X n (K n (s)))\ 
+\X n (s) - X n ( Kn (s))\\b(s,X( S )) - b n (s,X n ( Kn (s)))\ 
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and on further application of Young's inequality and A-l, this becomes 



(e„(s)A(s)> < (2C R + l)\e n (s)\ 2 + (2C R + l)\X n (s)-X n (K n (s))\ 2 

+-\b(s,X n (K n (s))) - b n {s,X n {K n (s)))\ 2 

+\X n (s) - X n (K n (s))\{c(\X(s)\ + \X(K n (s))\ + M + M n )}. 

Therefore raising power p/2 on both sides, we get 



\(e n (s),b n (sW < 2P- 3 {2(2C R + l)f \e n (s)\P + 2(2C R + l)*\X n (s) - X n ( Kn (s))\P 
+ \b(s,X n (K n (s))) - b n (s,X n (K n (s)))\ p 

+2\X n (s)-X n (K n (s))\2{c(\X(s)\ + \X(K n (s))\+M + M n )}*} 

for every s £ [to,u A v n Ft\- Hence one obtains the following estimate on using Holder's inequality, 



77 := (32T) E 2 1 E 



ruAv nR _ p 

'to 



< (32r)VJ 2(2/(i2) + l)§ jf E 



sup \e n {r AVnR)\ p lc R <f(R) 

to<r<s 



ds 



+2(2/(i?) + l)2E 



to 



+E 



\X n (s A f nJ? ) - X n (K n (s A v nR ))\ p ds 
\b(s A A v nR ))) ~ b n {s A v nR , X n (K n (s A ^w?)))| p <is 



+2P +l/2 Tc f 







E 


/' 

/to 



\X n (s A - X„(k„(s A v nR ))\Pds 



x v^TA 7 



(17) 



Now in order to estimate III of (fPoj) . one observes that for every i S {1,. . . ,m}, the assumption 
A-2 yields 

14(5)^ = |a i ( S ,X( S ))-4(s,X4/^( S )))P' < ^{cllenCs)!' + c||X n ( S ) - X„( Kn ( S ))| p 

+ |cT*( S , X n ( Kn ( S ))) - <7* ( S , X n ( Kn ( S )))| P } 
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which implies 



uAu nR 

Hc R <f(R)} I K(s)\ p ds 
to 



III := (4mT)V^E 
i=i 

< 3 p-1 (4mT) I! 7 2 |m/(l?)5 jf 

+m/(/2)*E 



E 



sup |e„(r A z/ nR )| p l {Cfl < /(jR)} 



/ |X„(s A i/ nJJ ) - X n (/c n (s A ^„i?))| p (is 
7t 



i=l 



K(s A ^ nfi ,X n (K„(s A fnfl))) - o-„(s A ^ nR ,X n (K n (s A v n R)))\ p ds 



to 



(18) 



Finally, to estimate IV, observes that for i £ {1, . . . , m}, an application of assumption A-2 yields 



(e n (s),4(s))|2 < |e n (s)|2|al(s)|2 = |e n (s)| 2 ^(s, - cr^s, X n ( Kn (s)))\ 2 



which by using Young's inequality becomes 

\(e n (s),a n (s))\^ < 3^{c|| en ( S )|f + c||e n (a)|§|X n (s)-X n (K n (a))|§ 

+|e n («)|§|<7 i (a,X n (K n (a))) - 4(s,X„( Kn ( S )))|l} 

'3C| + l) |e n (s)| p + c||X n (s) - X n ( Kn ( S ))| p 

+|<7 < (a, *„(«„(*))) - a n (s,X n (K n (s)))\p]. 



< 



3 2 



Thus 



3p — 4 p— 2 p — 4 ^ — -\ 

IV := 2 — m—T—c22 E 

i=i 



ruAu nR 

1 {C R <f(R)} / \{en(s),a t n (s))\^ds 
J t 



3(p-2) p-2 p-2 p-4 _ , , 

< 2-^-Z—m—T— c{m(3f(R)* +1 / E 



to 



sup |e n (r A v nR )\ p t {CR < f{R)} 

to<r<s 



ds 



+mf{R)m 



[ \X n (s A v nR ) - X n (K n (s A v nR ))\ p ds 
Jt 



i=l 



[ \a\s A ^ nR ,X„(K n (s A VnR.))) ~ <? n ( s A V nR ,X n {Kn{s A ^i?)))^^ 
7 to 



■(19) 
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Hence using (fT7|) . (JTHJ) and (fl"9j) in (fT5j) and then after applying Gronwall's inequality, one obtains 



E 



SUp \e n (t /WnR)\ P t{C R <f(R)} 
t <t<h 



<e^ T \2P- 2 E[\e n (t )\P) 



+K 2 E 

m 
i=l 



\X n (s f\ V nR ) - X n (K n (s A V n R))\ P ds 



\b{s A v nR ,X n (n n (s A v nR ))) - b n (s A v nRl X n (K n (s A v nR )))\ p ds 



to 



K(s A u nR ,X n (K n (s A v nR ))) - a l n (s A v nR ,X n (n n (s A ^„ij)))| p (is 



E 



\X n (s A i/ nii ) - X n (n n {s A f nR ))|Pds 



to 



(20) 



where constants K\ , Ki , K3 , -^4 and ET5 are appropriately defined and depend explicitly on f(R) 
but not on R and they are independent of n. Thus one can choose r/ sufficiently small and R 
sufficiently large such that for e > (however small), 



^o9^/ e q-p2K e 

— l q K < -, 77 r— — - < - and 

q 3' qrjP/iv-p) BP 3 



9 ,, P M Cr > f(R)) < % 

p/(q-p) v n J \ II 3 



Therefore, substituting equations (fT4"|) and (f20|) in (fT3j) and then using assumption A-3, A-4 and 
LemmaEl one obtains (|12p for p > 4. The application of Holder's inequality completes the proof. □ 



4 Proof of Main Result 

In this section, we shall prove Theorem [1] by considering SDDE ([TJ) as a special case of an SDE 
with random coefficients. In particular, we set 

b{t, x) := P(t, Y{t),x),a(t, x) := a(t, Y{t),x), b n (t, x) := /3(t, Y n (t),x),a n {t, x) := a(t, Y n {t),x) (21) 
for any t G [0, T] and x G M. d . 

Proof of theoremUl We shall prove the result using an induction method so as to show that, 

= (22) 



lim E sup \X(t) -X n (t)\ p 

«->oo (i-l)r<t<ir 
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for every i G {1, . . . , ./V} and for all p > 0. 

Case: t G [0,r]. When t G [0,r], then SDDE Q and Euler scheme (|2J) become SDE (jU) and 
scheme ([5]) with to = 0, t\ = r, X(0) = X n (0) = £(0) and coefficients given by ([2T]). Also one 
observes that A-l to A-5 hold due to C-l to C-3. In particular, A-l is a consequence of C-l with 
M = M n = 1 + V 1 , where * := sup te r i |(£(£i(*))> • • • , £(<5*;(*)))|> which is uniformly bounded 
due to the fact that £ G Cj}r . Further, A-2 is a consequence of C-2 since one observes that Cr 
is deterministic and in fact Cr = f(R), where f(R) = -K#l{o<.R<*} + ^Cr1-TR>#} an d therefore 
P(C R > f(R)) = for any R > 0. Finally, A-3, A-4 and A-5 hold trivially. Therefore, for i = 1, 
equation (|22p holds due to Theorem [2] and Lemma [1] for to = and t\ = r. 

For inductive arguments, we assume that when % = r, i.e. i G \{r — l)r, rr] for some r G 
{1, ... ,N — 1}, equation (|22p is satisfied and Lemma [T] holds for t$ = (r — l)r and t\ = rr. Then 
we claim that when i = r + 1, i.e. t G [rr, (r + l)r], equation ([22]) holds and Lemma [1] is true for 
to = rr and t\ = (r + l)r. 

Case: t G [rr, (r + l)r]. When t G [rr, (r + l)r], then SDDE ([I]) and its Euler scheme ([2]) become 
SDE flU) and scheme © with t = rr, ti = (r + l)r, X(t ) = X(rr), X n (t ) = X n {rr) and 
coefficients given by (|2~TT) . 

Verify A-l. Consider for t G [rr, (r + l)r] and x G M d , then due to C-l, 

|6(t,x)| + \a(t,x)\ = \/3(t,Y{t),x)\ + \a(t,Y(t),x)\ < G{M+\x\) 



and 



\b n (t,x)\ + |<7„(t,x)| = \/3(t,Y n (t),x)\ + |a(t,y n (t),x)| < G(M n + \x\) 

where M := 1 + sup rr<t <( r+1 ) r |l"(t)|' and M n := 1 + sup rT<4< ( r+1 ) T |Y^(t)| which are bounded in 

C p for every p > because of Lemma Q] and the inductive assumptions. 

Verify A-2. For every R > 0, |z| < and t G [rr, (r + l)r], from C-2, one obtains, 

(x - z, b(t, x) - b(t, z)) = (x-z, f3(t,Y(t),x) - P(t, Y(t),z)) < C R \x - z\ 2 
\a(t, x) - cr(t, z)\ 2 = \a(t, Y(t),x) - a(t, Y(t), z)\ 2 < C R \x - z\ 2 
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where the random variable Cr is given by 

oo 

C R := K r 1q r + ^2 K j+1 l {n . +1 \ Q . } 
j=R 

with Qj := {oj £ 0, : sup te r rT) ( r+1 ) T i \Y(t)\ < j}. One observes that Cr is an J> T -measurable random 
variable. Further one takes f(R) := Kr for every R > so that 

¥(C R > f(R)) = ¥(Cr>Kr) = 1-¥(Cr<K r ) 

< l-F(Sl R ) = w( sup \Y(t)\ > r) -> as R -> oo. (23) 

Verify AS. For R > and |x| < R, one observes that \Y(t) - Y n (t)\ 4 as E[|Y(t) - Y n (i)| p ] -»• 
for every p > when n — > oo. Therefore from C-3, 



sup |/3(t,y n (t),x)-/3(t,y(t),x)| p 4 and sup \a{t, Y n (t), x) - a(t, Y(t), x)\ p 4 0. 
\x\<R \x\<R 



Furthermore, one observes that sequences 



sup \/3(t,Y n (t),x)-P(t,Y(t),x)\ p \ and\ sup \a(t, Y n (t), x) - a(t, Y(t), x)\ p 

\x\<R J n >l I \x\<R 



n>l 



are uniformly integrable since they are bounded in C q for any q > 1, 



E 



sup \/3(t,Y n (t),x)-p(t,Y(t),x)r 
L |x|<R 



sup |/3(t,y n (t),x)-/3(t,y(t),x)r 

L|a;|<R 



< 2 P9_1 E 



E 

+ 2 P9-1 E 

< 6 p9-1gp<?|i +E[|y n (t)| Ip9 ] + \R\ pq } + 6 pq - 1 G pq {l + E[\Y(t)\ lpq ] + 

< & ,q ~ 1 G pq 2{l + K+ \R\ pq } 



sup |/3(t,y n (t),x)r 

|x|<R 



sup \/3(t,Y(t),x)\ pq 
\x\<R 



and similarly, for the sequence {swp\ x \ <R \a(t,Y n (t),x) — a(t, Y(t), x)| p } n >i- Therefore, 



E 



sup \\P(t, Y n (t),x) - f3(t, Y{t),x)\ p + \a(t, Y n (t),x) - a(t, Y(t), x)\p\ 



L M<-R 
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as n — > oo due to Dominated Convergence Theorem which also yields 



lim E 

n— >oo 



(r+l)r 



sup {\b(s,x) - b n (s,x)\ p + \a(s,x) - a n (s,x)\ p }ds 

\x\<R 



Thus A-3 is satisfied. 

Verify A- 4- This holds due to the inductive assumptions. 
Verify A- 5. This follows due to the inductive assumptions. 

Finally Lemma[Uholds for to = rr and t\ = (r+l)r due to the fact that A-l holds for t G [rr, (r+l)r] 
and x G M. d and A-5 is true for to = rr - This completes the proof. □ 



5 Rate of Convergence 

In this section, we shall recover the rates of convergence of the Euler scheme j2j) under different set 

of assumptions. First we state below the relevant assumptions. 

C-4. There exist constants C > and l\ > such that for any t G [0, T], 

(x - x ,/3(t,y,x) - f3(t,y,x'))V \a(t,y,x) - a(t,y,x')\ 2 < C\x-x'\ 2 

\P(t,y,x) - P(t,y ,x)\ 2 + \a(t,y,x) - a(t,y ,x)\ 2 < C{1 + \y\ h + \y \ h )\y - y'\ 2 

for all x,x' G R d and y,y' G R dxk . 

C-5. There exist constants C > and I2 > such that for any t G [0,T], 
\P(t,y,x) - (3(t,y,x')\ 2 < C{1 + \x\ h + \x'\ h )\x - x'\ 2 
for all x, x G R d and y G M. dxk . 

In Corollary [H we prove that under C-l and C-4, the rate of convergence of Euler scheme ([2]) 
is one-fourth. On the other hand, in Corollary [2] we prove that if one makes assumption C-5 in 
addition to C-l and C-4, then the classical rate of convergence (one-half) can be recovered (see also 
d, EH, EH], M and [20]). 

Corollary 1. Suppose C-l and C-4 hold, then for all p > 0, the Euler scheme ([2]) converges to the 
true solution dU in CP -sense with rate 1/4, i.e. 
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E 



sup \X(t)-X n (t)\ p 

0<t<T 



< Cn-i 



where the constant C > does not depend on n. 



Proof. In equation (|14p . we choose rj = n 2 , R = n 2 ^-^ , q > 2p > 4, then 



q _£ q-p 2K 2{q-p)K _ £ 
n 2 , 77 r-^— = n 2 and 



P 



•p/(q-p) Rp 



pI (q-p) ' 



\C R > f(R)) = (24) 



where we choose Cr = C and f(R) = C. Then one can show that, 



E 



sup \X(t) - X n (t)\ p 

(k-l)T<t<kr 



< Cn~4 



(25) 



for each k G {1, . . . ,N}. As seen before, SDDE (UJ corresponds to the ordinary SDE @ with 
X(0) = X n (0) = ^(0), Y(t) = Y n {t) = £{5{t)) for all t G [0,r]. Therefore, one obtains the rate as in 
(|25p for k = 1, which is 1/4 instead of 1/2, due to the last term of (|17p . that is also the last term 
in (|20p . the estimate in (jlip and the fact that 



E 



\b{s A u nR ,X n {K n (s A v nR ))) - b n (s A f nJ? ,X n (K n (s A v n R)))\ p ds 







and 



E 



/' 

Jo 



\o- l {s A i/ nfi ,X n (K n (s A fnR))) - o- l n {s A v nR ,X n (K n (s A twj)))| p cfe 



= 



for every i G {1, . . . , m}. 

One then assumes that (]25p holds for k = r, i.e. t G [(r — 1)t, rr] for some r G {1, . . . , N — 1}, 
so as to show that it also holds for t G [rr, (r + l)r]. Due to inductive assumption, the initial data 
of SDE (01) satisfy 



E 



\X(rr)-X n (rr)\P 



(26) 
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and due to (|2ip . C-4, Lemma [T] and inductive assumptions, one obtains the following estimates, 



E 



E 



(r+l)r 



\b(s A v nR ,X n (K n (s A v nR ))) - b n (s A v nR ,X n (K n (s A v nR )))\ p ds 



(r+l)r 



< 



(r+l)r 



< 3 £ 2 i v / ( 1 + 2 ^) 



|/3( s ,y( s ),x„( Kn ( s ))) - /3( s ,y n ( s ),x n (K n ( s )))|fi {s <^ } ^ 



e[(|i + ly^)!' 1 + |y n ( s )|' 2 )f|y( s ) - Y n ( s )\ p ]ds 

(r+l)r 



E[|Y(a)-y n (8)|2p]dfl 



< Cn _ 5. 



(27) 



One similarly estimates 



E 



(r+l)r 



A A f n fl))) - °"n( s A u nR, X n (K n (s A l^i?))) I^S 



< Cn~2 (28) 



for every ie{l,...,m}. Therefore, ([25]) holds due to ([20]) by taking into account (pi]), ([Ml), ([26]) . 
([271) and ([28]) for i G [rr, (r + l)r] and p > 4. Then the application of Holder's inequality completes 
the proof. □ 

Corollary 2. Suppose C-l, C-4 and C-5 hold, then for all p > 0, the Euler scheme ([2]) converges 
to the true solution ([!]) in C p -sense with rate 1/2, i.e. 



E 



sup \X(t)-X n (t)\ p 

0<t<T 



< Cn~2 



where the constant C > does not depend on n. 

Proof. In order to prove this, one observes that the estimates of the rate of convergence of terms 
of // in (117p is improved if C-5 is used along with C-4. For this, one could replace the right hand 
side of ([T6]) by the following, 



(e n (s),6 n (s)) = (e n (s),b(s,X(s)) - b(s,X n (s))) + (e n (s), b(s, X n (s)) - b(s,X n (K n (s)))) 

+ (e n (s),b(s,X n (n n (s))) - b n (s,X n (K n (s)))} 
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which on using (|2ip becomes 



(e n (s), b n (s)) = (e n (s),P(s, Y(s),X(s)) - f3(s, Y(s),X n (s))) 

+ (e n (s),P{s, Y(s),X n (s)) - p[a, Y(s),X n (K n {s)))) 

+{e n (s),j3(s, Y(s),X n ( Kn (s))) - (3(s, Y n (s) , X n (K n (s)))) . 

By using C-4 and Cauchy-Schwartz and Young's inequalities, one obtains 

(e n (s),b n (s)} < (C + l)\e n (s)\ 2 + l\f3(s,Y(s),X n (s)) - (3(s,Y(s),X n (K n (s)))\ 2 



+-\(3(s,Y(s),X n (K n (s))) - l3(s,Y n (s),X n (K n (s) 



which due to C-4 and C-5 yields 



(enOOAOO) < (C + 1)M*)| 2 + ^{1 + \X n (s)\ h + \X n ( Kn (s)))\ h }\X n (s) - X n ( Kn (s 

+^{1 + \Y(a)\^ + |y n (*)|'i}|y00 - Y n (s)\ 2 . 



Therefore, one estimates II in (|17p by the following, 



II := (32T) £ 2 1 E 



i-uAu„ R _ 
'to 



< (96T)V^ (C + l) 



P-2 



— / E 

rr 



sup |e n (sj 



ds 



rr 

p 



{1 + |X n ( S )| i2 + |X n ( K „( S )))| i2 }l|X„( S ) - X n (K n (s)))\fds 
{1 + ^(s)!' 1 + \Y n (s)\ h } P i\Y(s) - Y n (s)\Pds 



which on the application of Holder's inequality and Lemma [T] yields 



II < i'HiVyV' \ (C7 + 1)V / IE 



{ ( _y h r +l,T 



sup |e n (s)| p 

rr<t<s 



ds 



E 



|X n ( S )-X n ( Kn (s)))| 2 P +WE \Y(s)-Y n (s)\ 2 P 



ds 
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where K is a constant which does not depend on n. Therefore, in (|15p . one replaces the estimate 
(|17p of II by the estimate ()29f) whereas I// and /V remain the same. Then the first term of (|29p is 
incorporated with similar terms of III and IV so as to apply Gronwall's inequality. Furthermore, 
one observes that the second and third terms of (|29p are of (improved) order 0(n -p / 2 ) due to 
Lemma [2] and inductive assumptions respectively. Thus, by adopting same arguments as adopted 
in the proof of Corollary [1] with 7] = n~2, one obtains the desired rate (1/2). □ 

6 Numerical Examples 

In this section, we shall illustrate our findings with the help of numerical examples. We consider 
the following SDDE given by 



dZ(t) = [aZ{t) + b{Z(t - r)} h ]dt + [/3t + f3 2 Z(t) + (3 3 {Z(t - r)} h ]dW(t) for t € [0, 2r] (30) 

with initial data £(t) when t G [— r, 0], r > is a fixed delay, > 0, a, 6, ^1,^2,^3 G R. 
When i G [0, r], then explicit solution of SDE (f30|) is given by, 

Z(t) = $ ,t{m + j *oi{b{t(* - T )} h - Wi + /? 3 {e(s - r)} h )}ds 

+ J *oiifo + - r)} l *}dW(s)}, (31) 

where, 

$o,t = exp ({a - + . 

And when i G [r, 2r] , then the explicit solution is given by 

Z(t) = * T)t |z(r) + £ " ^ " Wi + - T)} l *)]ds 

+ jf + h{Z{s - r)} l *)dW(s) 1 (32) 
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where, 



* T ,t = exp ({a - ^}(t - r) + ft(W(t) - W(r 



For Euler scheme, we divide the interval [0, r] into sub-intervals of width for some positive 
integer N. Then solution at kh-th grid is given by 



Z n {(k + l)h) = [aZ n {kh)+b{i{kh-r)} h ]h 

+ [ft + /3 2 Z n (M) + /3 3 {^(^ - r)}' 2 ] + - W{kh)) 

with £ (i) for t G [-r, 0], fc = 1, . . . , 2^. We use this solution to find the solution of the Euler scheme 
in the interval [r, 2r] . 

We take values of the parameters of SDDE (|30|1 given in Table [TJ From Table El we observe 
that as step size is decreased, the error is decreased too. In Figure [TJ the reference line has a slope 
of -0.5. 

Table 1: Values of Parameters 



V 


r 


a 


b 


ft 






£(t) for t G [-r, 0] 


2 


1 


-8 


4 





1 


1 


t+1 



Table 2: Effect of decreasing step size on Error 





VE(\Z(T) - Z n (T)\ 2 ) 


h 


h = 1/2, Z 2 = 1/2 


h = l,h = l 


h=2,l 2 = 3 


2-9 


0.0002332678715832590 


0.0001857299504105190 


0.0002812168347124360 


2-io 


0.0001171164266162860 


0.0000911946799085898 


0.0001308779599825470 




0.0000605972526234064 


0.0000450951550271487 


0.0000652194083354262 


2-12 


0.0000314672425957171 


0.0000223376219662284 


0.0000321943295689175 


2 -13 


0.0000159827178333670 


0.0000111954499489959 


0.0000156672183324178 


2-14 


0.0000079146081811971 


0.0000052537027007630 


0.0000057778221373887 


slope 




-0.51157 


-0.54615 
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Figure 1: Strong order of convergence of the Euler scheme. 
(a)Zi = l,Z 2 = l (b)Zi=2,Z 2 = 3 
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